Non-local Coulomb correlations in metals close to a charge order insulator transition 
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The charge ordering transition induced by the nearest-neighbor Coulomb repulsion, V, in the 
1/4- filled extended Hubbard model is investigated using Cellular Dynamical Mean-Field Theory. 
We find a transition to a strongly renormalized charge ordered Fermi liquid at Vco and a metal-to- 
insulator transition at Vmi > Vco- Short range antiferromagnetism occurs concomitanly with the 
CO transition. Approaching the charge ordered insulator, V < Vmi, the Fermi surface deforms and 
the scattering rate of electrons develops momentum dependence on the Fermi surface. 
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Strongly correlated electronic systems in low dimen- 
sions exhibit novel quantum states of matter with un- 
conventional electronic properties. Among the fascinat- 
ing phenomena observed in these systems is the Mott 
metal-insulator transition induced by the local Coulomb 
repulsion energy [l| which occurs, for instance, in transi- 
tion metal oxides and layered organic superconductors 
[3, 0]. Important progress in the description of the 
Mott transition [1] has been achieved with the develop- 
ment of dynamical mean-field theory (DMFT)f5|, a local 
theory which predicts that close to the Mott insulator 
(MI) , a quasi-particle peak and Hubbard bands coexist in 
agreement with photoemission spectra of transition metal 
oxidesQ and optical spectra of k-(BEDT-TTF)2X 0. 
Non-local correlations can modify this scenario opening 
a pseudogap in the spectral density as Cellular Dynami- 
cal Mean-Field Theory (CDMFT) Q and the Dynamical 
Cluster Approximation (DCA) [9| have found consistent 
with the pseudogap phase of high-Tc superconductors. 

The 1/4-filled family of layered organic molecular 
superconductors of the a, /3" and 0^BEDT-TTF)2 
X types display rich phase diagrams [lOj and anoma- 
lous metallic behavior 11] such as the Drudc peak ab- 
sence in the optical conductivity of ^"-(BEDT-TTF)2 
SF5 CH2 CF2 SO3 [l3| ■ In these quasi- two-dimensional sys- 
tems both the on-site Coulomb repulsion, U, and the 
nearest-neighbor Coulomb repulsion V are substantial. 
In spite of its importance, the electronic properties of 
metals close to a Coulomb-driven charge ordered insu- 
lator (COI) transition remain poorly understood. Tak- 
ing the extended Hubbard model as the minimal strongly 
correlated model which includes the competition between 
itinerancy and localization due to charge ordering we 
study the properties of metals close to a COI. An impor- 
tant theoretical challenge is to understand the evolution 
of the spin degrees of freedom across the charge ordering 
transition and, in particular, whether spin and charge 
order together or not. 

In this Letter, we describe how the COI transition 
driven by the Coulomb repulsion, V, in the extended 



Hubbard model occurs based on the CDMFT approach, 
which is a strong coupling method which includes both 
on-site and short range non-local correlations. We find 
a transition to a strongly renormalized charge ordered 
metal (COM) at Vco and a metal-insulator transition 
at Vmi- We find that the spectral density is suppressed 
on some regions of the Fermi surface due to the momen- 
tum dependence developed in the electron scattering rate 
close to the COI. Concomitanly with the CO transition 
an antiferromagnetic (AF) interaction between spins at 
charge rich sites of the COM is dynamically generated by 
the Coulomb interaction, V. 

The simplest strongly correlated model that captures 
the competition between charge ordering and the homo- 
geneous metal is an extended Hubbard model on a square 
lattice at 1/4-filling: 



(1) 



where a is the spin index and and Cia denote the 
electron operators, t is the transfer integrals between 
nearest- neighbor sites of the square lattice. The Coulomb 
interaction is such that: Vij = V \i i and j are nearest 
neighbor sites and Vij = otherwise. The system with 
no interactions {U = V = 0) has a dispersion with band- 
width W = 8i. We fix C/ = IQt and analyze the CO 
transition due to V at 1/4-filling. 

In the CDMFT method, the infinite lattice is divided 
into identical clusters of size Nc = ^ which are treated 
exactly as shown in Fig. [T] (a). The Coulomb interaction 
between electrons in different clusters is treated at the 
mean-field level. The quantum cluster is embedded in 
a non-interacting bath which is solved self-consistently. 
Results are presented for paramagnetic phases unless ex- 
plictly stated. The effective Anderson model reads 
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where tij = —t denote hopping matrix elements between 
neighbor sites inside the cluster with i and j running from 
1 to Nc. Following Ref. bath levels are grouped into 
a subsets. Correspondingly the energies, ej?, refer to the 
k energy level of the a subset. The bath-cluster couplings 
between a bath level k of bath a and the cluster site i 
are: FJ?-. In practice, we truncate the bath to A^b = 8 
sites which are grouped into two subsets: a = 1, 2 with 
fc = 1, .., 4 states each. We take iVc = 4 cluster sites. Lat- 
tice symmetries (allowing for checkerboard charge order- 
ing) reduce the bath parameters to 4 independent bath 
levels: ej? — e5?_|_2) and hybridization bath parameters: 
F^^ = 4iFf , with Ff = Ff+2- Thus, we use 4 indepen- 
dent bath levels and 4 independent hybridization param- 
eters to fully describe the bath. After extracting the self- 
energy of the cluster, Sy (iaj„), for an arbitrary choice 
of tiS and F^i's the lattice Greens function is found by 
imposing the superlattice periodicity: 

GijiiUn) = y^((i^,t+M-y^ Vik{nk))dij-tij{K)-Y.ij{iujn)) 
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where tij (K) is the Fourier transform of the cluster hop- 
ping amplitudes and K the momentum wavevector in 
the reduced Brillouin zone of the superlattice. Mat- 
subara frequencies, ujn = (2n + 1)tt/P, with an inverse 
temperature: (3 — 200/i are used, (n^) is the average 
electronic occupation in nearest-neighbor sites in neigh- 
boring clusters. The one-dimensional version of model 



(H]) has been previously studied with CDMFT[1J|. The 
two-dimensional Hubbard model on a frustrated square 
lattice close to half-filling (with small degree of frus- 
tration), leads to an intricate competition between AF, 
Mott insulator, unconventional superconductivity and 
metallicitv[l5|. 

The resulting phase diagram of the extended Hub- 
bard model based on CDMFT is shown in Fig. [1] As 
V is increased the system undergoes a sharp transition 
from an homogeneous metal (M) to a strongly renor- 
malized charge ordered metal (COM) at Vco ~ l-3i- 
The jump in the order parameter of the CO transition, 
ni2 = {ni) — (712) is evident from Fig. [T] (b). Initially, 
for V < Vco, the system is weakly correlated as the 
quasiparticle weight Zi ~ {1 — 9ImI]ii(iu;„)/9a;„|i^=o) in 
charge rich (Zi) and charge poor {Z2) cluster sites is 
Zi ~ Z2 = 0.7 consistent with slave-boson theory [Tij]. 
The weakly correlated metal becomes a strongly renor- 
malized Fermi liquid once charge ordering sets in with 
Zi — 0.1 at V = 1.5t indicating localization of electrons 
at every other site of the lattice due to CO. At Vmi ~ 2i a 
gap opens in the spectral DOS indicating a MI transtion. 
Previous direct Lanczos calculations on 20-site clusters 
predict a MI transition at Vmi ~ 2.2t and indicated a 



FIG. 1: (Color online) Ground state phase diagram of the 
1 / 4-filled extended Hubbard model on the square lattice from 
CDMFT calculations. In (a) the tiling of the two-dimensional 
lattice used in CDMFT is shown whereas in (b) the V- 
dependence of site occupations, Ui, is plotted for U = lOt. 



COM for V > Vco ~ l-5i in remarkable agreement with 
_<^UT present results despite the small size of the cluster 
used within CDMFT. 

Dynamical and short range correlation effects are ana- 
lyzed based on cluster quantities: Gijiiun) and Y^ijiiujn). 
Figs. [5] and [3] show how diagonal and off-diagonal self- 
energies are enhanced for V > Vco being most significant 
between charge rich sites (see Si3(ia;„) in Fig. [3]). Con- 
sistent with the self-energy enhancement, Green's func- 
tions are suppressed at finite LUn for V > Vco- For 

V > 2t, ImGii(iw„ — > 0) oc — w„, displaying insulating 
behavior [l3| • This contrasts to the metallic- like behavior 
of ImG'ii(iw„) for V < Vmi indicating a MI transition at 
Vmi = 2t. Fig. [3] shows ImSy (zci;„ — * 0) indicating 
that the opening of the gap is due to the large ReEi3(?0"'") 
developing close to the COI (see Fig. [3]). This interac- 
tion induced by Coulomb interaction produces a band 
sphtting at the Fermi energy^. Note that as LRO-AF 
is not allowed, the opening of the gap at Vmi can only 
be due to short range correlations. We also find that at 

V w 2.5t, ImEii(i0+) and ImEi3(i0+) diverge. 

The spectral DOS, Ai{uj), at site i is shown in Fig. 
m indicating metallic behavior for V < Vmi while a 
gap is apparent at Vmi ~ 2i. The metallic phase can 
be further explored through the momentum resolved 
spectral DOS, ^(k, w) with k being the momentum 
in the 1st Brillouin zone of the real lattice. Lattice 
self-energies are extracted by periodizing cluster self- 
energies: S(k,tj) = ■^^Y.i'^iii'^) + ^i2{c^)cos{k^) + 
£14(0;) cos(fcj,) -I- £24(^)003(^2: -I- ky)/2 + Y,i^{ijj)cos{kx — 
ky)l2^ with distances taken in units of the lattice pa- 
rameter. The evolution of A{k, 0+) is displayed in Fig 
[S] which compares V = t {a, metal with no CO) with 

V = IM (a COM close to the COI). For V ^t, A(k, 0+) 
describes a somewhat renormalized Fermi liquid and the 



3 




FIG. 2: (Color online) Enhancement of on-site self-energies 
across the CO and MI transitions with V. In (a) Eii(iLj„) 
and Gii{iuj„) which correspond to the charge-rich sites are 
plotted and in (b) the same for the charge-poor sites. 
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FIG. 4: (Color online) Evolution of the DOS with the 
Coulomb repulsion, V, across the CO transition. Ai{ijj) is 
plotted for different V showing a gap opening at the COI 
transition Vmi ~ The vertical line denotes the position of 
the chemical potential. 




FIG. 3: (Color online) Dynamically generated off-site interac- 
tion in the COM. Non-diagonal self-energies are plotted with 
V. Self-energies T.ij{iuj„) between two next-nearest neighbor 
charge rich sites are plotted in (a) while between nearest- 
neighbor sites are displayed in (b). 



non-interacting Fermi surface shape. For V ^ Vmi the 
Fermi surface is deformed [l^ and the scattering be- 
comes momentum dependent. 'Cold' regions develop 
around momentum k w (±0.4, =F0.4)7r whereas 'hot' re- 
gions occur around k « (±0.8, ±0.8)7r. This arises from 
the momentum dependence of the electron scattering 
rate: 1/rk = — 2ImI](k, 0"^) approaching the COI. In 
the top part of Fig. [5] the decrease of l/rk along the 



(0, 0) (tt, — tt) direction is sharper for larger V while 
along (0,0) — !■ (tt, tt), 1/rk remains large. The momen- 
tum dependence of 1 /rk leads to the anisotropic momen- 
tum dependence of ^(k, 0"^). Although the precise posi- 
tion of 'cold' and hot 'regions' can vary for larger Nc (as 
the momentum dependence of the self-energy is restricted 
by the symmetries of the small Nc = 4 cluster), the mo- 
mentum dependence of 1 /rk and Fermi surface deforma- 
tion are robust features derived from CDMFT close to the 
COI. This is in contrast to single-site DMFT|20|] which 
predicts a direct transition from a metal to a strongly 
correlated Fermi liquid with isotropic scattering rate for 
any V. 

Magnetic properties are analyzed by computing the 
static spin correlation function, (SizSjz), between dif- 
ferent sites in the cluster. The spin interaction between 
charge-rich sites, (SizSsz), is AF for all V as shown in 
Fig. [6l At the charge ordering transition there is a sud- 
den enhancement of (SizS^z)- Part of it is attributed 
to the charge transfer which increases the magnitude of 
spin moments at the charge rich sites. An AF interaction 
is expected in the strong coupling limit, U/t, V/t — > oo, 
which lead to an effective J = Af^ /9V^ due to fourth or- 
der 'ring' exchange processes 3 21 1. The possibility of 
having long range order antiferromagnetism (LRO-AF) 
induced by this AF interaction is explored by allowing 
AF solutions within CDMFT. As shown in Fig. [H the 
staggered magnetization: m = 1/A^c X]i(~l)'("'iT ~ 
becomes finite at the CO transition, therefore LRO-AF 
develops concomitanly with CO. This does not discard 
the possibility of having a paramagnetic COM at finite 
T as shown previously [2^]. This is because the T-scale as- 
sociated with charge ordering Tqo is governed by V while 
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FIG. 5: (Color online) Momentum dependence of spectra 
close to the MI transition driven by V. Top (bottom) fig- 
ures show the evolution of 1/rk (A(k, 0^)) with V. Left plots 
show results for V = t while right plots are for V = 1.8t 
(U = lot). Regions with large 1/rk and suppressed j4(k, 0"'') 
are found close to the COI {V — 1.8t). Blue (red) color corre- 
spond to (zero) largest intensities being for bottom left plot: 
(blue), 1.8 (yellow) and 2.6 (red) and for bottom right: 
(blue), 1.5 (yellow) and 2.1 (red). 
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FIG. 6: (Color online) Antiferromagnetic interaction for 
U = lot. In (a) a sharp increase of the spin-spin correla- 
tion, (SizSjz) between charge rich sites occurs at Vco ~ l-St. 
In (b) the staggered magnetization shows that the onset of 
LRO-AF coincides with CO. 



AF ordering is controlled by Taf ^ J << V. Hence, the 
CO paramagnetic metal with a momentum dependent 
l/rfe should exist in the T-range: Taf < T < Tqo- 

In summary, as the COI transition driven by V is ap- 
proached the system evolves from a weakly renormal- 
ized Fermi liquid to a metal in which the scattering rate 
becomes momentum dependent which suppresses spec- 



tral weight on some parts of the Fermi surface. Short 
range AF is dynamically generated by V at the CO tran- 
sition. Angular resolved photoemission experiments on 
1 /4- filled 9 and /3" layered organic molecular compounds 
are highly desirable to test CDMFT predictions. 
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